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Abstract 

Based on the first-principles calculations, we predict a novel half-metallic surface state in layered 
bulk insulator NaCo02, with tunable surface hole concentration. The half- metallic surface has a 
single fermi surface with a helical spin texture, similar to the surface state of topological insulators, 
but with the key difference of time reversal symmetry breaking in the present case. We propose 
realization of topological superconductivity and Majorana fermions when the half-metallic surface 
states are in proximity contact with a conventional superconductor. 
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The possibility of realizing Majorana fermions in condensed matters as emergent excita- 
tions has stimulated great current interest. Majorana fermions are particles which are their 
own antiparticles They constitute only half of a usual fermion, and obey the non-Abelian 
statistics |2j, which is the key ingredient for the fault-tolerant topological quantum computa- 
tion jjjl. The Majorana states associated with zero energy mode have been predicted to exist 
in various systems , nevertheless, their experimental realizations remain challenging 

due to the requirements of extreme conditions such as strong magnetic field, low temper- 
ature, or ultra-clean samples. In this letter, we will show that NaCoC^/superconductor 
hetorestructure is a new and simple platform for realizing such exotic states. 

The simplest Majorana bound state is associated with a vortex core or edge in a two- 

nn 

dimensional (2D) topological p x + ip y superconductor [5|, [6|, which has a full pairing gap in 
the bulk but with gapless chiral edge states (which consists of Majorana fermions). Unfortu- 
nately, such topological superconductors are very rare in nature, this leads to the proposals 
of possible "induced" p x + ip y order parameter through the proximity effect, particularly 



for materials in contact with the simplest s-wave superconductors [iJl^ L-The proposal by 
Fu and Kane |7j is to use the surface state of topological insulators 13l-ll5|. where the spin 
degeneracy is removed and yet strong proximity effect can be expected due to the charac- 
teristic helical spin texture of Dirac-type surface state. The experimental setting can be in 
principles obtained from laboratory, however, since most known three-dimentional topolog- 
ical insulators 15Ml8| up to now are not good bulk insulators and important surface states 
may overlap with bulk states, experiments have to wait for development of well-controlled 
clean samples. It was proposed recently that semiconductor quantum wells with Rashba 
type SOC in proximity to s-wave superconductor will produce similar effect [a-llll]. This 
may lower the experimental threshold, since well-controlled samples are available nowadays. 
Both proposals are encouraging, while experimental obstacles still remain. First, and most 
importantly, magnetic insulating layers or strong external magnetic field are required to 
break the time reversal symmetry, which is not easy to implement experimentally; second, 
the Fermi surfaces in both cases are too small, and fine control of chemical potential is 
difficult for semiconductors in contact with a superconductor. 

To avoid the complication required by breaking time reversal symmetry (such as the 
usage of magnetic insulating layers or magnetic field), the most natural way is to start from 
a magnetic compound, rather than non-magnetic topological insulators or semiconductor 
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FIG. 1: (Color online) The structure of NaCo02 slab consisting of 5 unit layers. After structure 
optimization, the surface Co-0 bond-lengths are slightly modified (by 0.07A). 



quantum wells. On the other hand, however, the following ingredients have to be satisfied 
over a wide energy or doping regime, in order to induce the p x + ip y superconductivity 
through proximity effect: (1) two dimensional metal with a single (or an odd number of) 
Fermi surface; (2) strong enough helical spin texture arising from spin-orbit coupling (SOC). 
Those conditions all together suggest that a 2D single band half-metal materials with strong 



enough SOC will be the best candidate 12|, |l9j. Following this strategy, we propose in 



this work that NaCo02 is such a unique compound, which satisfies all those conditions 
simultaneously. 

Bulk Na x Co02 is a well known layered compound crystallizing in planar triangle lattice 
with Co atom coordinated by oxygen octahedrally (Fig.l). The Na atoms are interpolated 
between C0O2 layers, and its concentration x can be systematically tuned, resulting in 
complex magnetic and electronic phase diagram 

H. In p artiCula , the uncon— 
superconductivity around £=0.35 [2l|, and the layered antiferromagnetic (AF) phase around 



x > 0.65 



22M24| . have drawn lots of attentions. For our purpose, however, we consider the 



stoichiometry NaCo02 (i.e. x—1.0) and its surface state. 

NaCo02 single crystal (with R3m symmetry) is experimentally available 25], and it is 
a simple band insulator with ba nd g ap more than l.OeV {25)]. Its bulk insulating behavior 
has been indicated by transport 26 1 and NMR {2^ measurements. Electronically, the Co- 
3d states splits into ti g and e g manifolds under oxygen octahedron crystal field, and all t^ g 
states are fully occupied (by 6 electrons) in the case without Na deficiency (corresponding to 
nominal Co 3+ case), resulting in an band gap between t 2s and e g j^Q. The layered crystal 
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FIG. 2: (Color online) The calculated band structure and density of state (DOS) for NaCo02 slab 
with surface hole concentration y=0.3. (a) The band structure of non-magnetic (NM) state. The 
projections to the a± g and e' g states of surface Co are indicated as red circle and blue triangle, 
respectively, (b) and (c) The corresponding DOS for the surface and bulk Co sites, respectively. 
The character of insulating bulk with metallic surface is seen, (d) and (e) The spin-resolved band 
structure of ferromagnetic state with the V z defined as the exchange splitting of surface a\ g state. 

structure guarantees that samples can be easily cleaved (Fig.l), and two kinds of termina- 
tions, with or without top-most Na layer, may be realized. Nevertheless, since the top-most 
Na 1+ ions are highly mobile, its concentration can be tuned depending on experimental con- 
ditions, resulting in surface hole doping (indicated as y) but without modifying the surface 
structure (the C0O2 layer) significantly. In the extreme case, if all top-most Na are absent, 
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0.5 hole will be introduced (i.e. y=0.5). We study the (001) surface of NaCo02 by using the 
first-principles calculations based on the plane-wave ultra-soft pseudopotential method, and 
the generalized gradient approximation (GGA) for the exchange-correlation functional. A 
slab consisting of five NaCo02 unit layers thickness (Fig.l) and 20 A vacuum region is used 
for the surface study, and the SOC is included from the fully relativistic pseudopotential. 
The cutoff energy for wave-function expansion is 30 Ry, and we use 12x12 k-points mesh 
for the Brillouin zone sampling. The calculations are well converged with respect to above 
settings. The calculations are further supplemented by the LDA+Gutzwiller method [281 ]. 
in which the density functional theory is combined with the Gutzwiller variational approach 
such that the orbital fluctuation and kinetic renormalization coming from correlation are all 
self-consistently treated. 

By optimizing the surface (slab) structure without top-most Na layer, our calculations 
confirm that the surface C0O2 layer remain well defined with only slightly (about 0.07A) 
modification of surface Co-0 bond length (see Fig.l). Another key result obtained from this 
calculation is that the hole only goes into the surface C0O2 layer, while keeping all other 
C0O2 layers inside the bulk insulating. There is only one band cross the Fermi level (Fig.2), 
which comes from one of the t<i g states of the top-most Co sites, while all other Co-t2 9 states 
inside the bulk are fully occupied. It is therefore effectively a system with insulating bulk 



13j-[l5|. The key difference lies in the 



but metallic surface, similar to topological insulators 
fact that the surface state considered in this work breaks the time reversal symmetry, in 
contrast to the case of topological insulators. In reality, the surface hole doping can be 
tuned by modifying the surface Na concentration, by interface charge transfer, or simply 
by gating. In the case of proximity effect with a superconductor, such surface doping can 
naturally arise since many s-wave superconductors, like Al, Pb or Nb, are simple metals. In 
such case, the insulating NaCo02 bulk can even serve as substrate simultaneously. In the 
following discussions, we will therefore neglect the top-most Na atoms and use the virtual 
crystal approximation to simulate the surface doping effect. 

Fig. 3 shows the magnetic properties of the NaCo02 surface as function of hole concentra- 



tion y. The spin polarization is energetically favored for all hole concentration 
in close analogy to the layered AF phase of Na x Co0 2 bulk (x ~0.75) [20, 



(0 



22 



5< y <0), 



29], where 



each C0O2 layer (with ~0.25 hole) orders ferromagnetically and the in-plane ferromagnetism 
contributes mostly to the energy gain with relatively weak inter-layer AF coupling. The sta- 
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FIG. 3: (Color online) The calculated magnetic properties of NaCoO"2 surface as function of surface 
hole concentration y. (a) The stabilization energy of ferromagnetic (FM) state with respect to non- 
magnetic (NM) state, (b) The magnetic moment of surface C0O2 layer. 

bilization of ferromagnetism at NaCoC>2 surface can be intuitively understood from Stoner 



instability (also similar to the discussion addressed for the layered AF bulk Na .75CoO 2 |29|). 
Due to the elongation of oxygen octahedra around Co sites along c-axis, the Co-t2 g states will 
further split into a\ g and e' g manifolds, with a\ g higher in energy. The a\ g state has mostly 
3d^ z 2_ r 2 orbital character (with z defined along c), whose in-plane dispersion is relatively 



weak and "M" shaped as observed in ARPES experiment [30j . As shown in Fig.2(a), the 
wide region of nearly-flat-bands around the valence band maximum will produce a sharp 
DOS peak near the Fermi level (see Fig.2(b)), which leads to the Stoner instability, and 
favors the ferromagnetic (FM) ground state for the surface. This mechanism is further sup- 
ported by the electronic structure of FM solution (Fig. 2(c) and (d)), where only the ai g state 
is strongly spin-polarized and the polarization of e' g states is small. In the FM state, the 
spin-polarization is strong enough to make the top-most C0O2 layer a half-metal, as shown 
in Fig.2 and Fig.3, where the calculated magnetic moment exactly equals to the number of 
holes. There is only one spin-channel of a\ g state crossing the Fermi level, resulting in a 
half-metallic surface state with single sheet of Fermi surface. 

We have to be aware of the effect of electron correlation beyond the generalized gra- 
dient ap pro ximation (GGA) for the exchange- correlation potential. The LDA+Gutzwiller 



method 



has been shown to be a powerful tool to take into account the correlation ef- 



fect, and reproduce correctly the magnetic phase diagram of bulk Naa-CoO^ [29j. We have 
supplemented the LDA+Gutzwiller calculations for the surface, and find that the surface 
ferromagnetism is further stabilized (by about 8meV/surface Co for ?/=0.3). In fact, in the 
study for the bulk Na.j;Co02 (x >0.6), both GGA and LDA+Gutzwiller give qualitatively 



the same result, well compared to experiments [29[. The AF state of bulk Na x Co02 was 
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FIG. 4: (Color online) The Fermi surface and the helical spin texture of NaCo02 surface (for y=0.3) 
obtained from (a) first-principles calculations, and (b) evaluation of effective k-p model (Eq.l, see 
text part) with parameter ,4=0.28 eVA 2 , 5=0.53 eVA 6 , C=-2.4 eVA 6 , a=-0.066 eVA, V z =0.22eV. 
The in-plane components of spin are indicated as arrowed lines, while the perpendicular component 
(pointing to out-of-plane) is not shown. 

observed for x > 0.65 with maximum T c around 25K for x ~0.8, the layered ordering; with 
spin orientation perpendicular to the plane was confirmed by neutron experiments 23|, |24| , 
and its half-metallicity of in-plane electronic structure is also supported by the ARPES mea- 
surements 31, 32|. Considering the layered nature and the similarity between the bulk AF 
phase and the surface, we conclude that a single band half metal can be realized at NaCo02 
surface with tunable hole concentration. 

Turning on the the SOC, the up and down spin bands will couple, while the characteristic 
single sheet of Fermi surface still remains. In the presence of a surface, the asymmetrical 
surface potential will produce the Rashba type SOC (which is automatically included in 
the first-principles calculations). It turns out that the Rashba SOC plays important roles, 
resulting in the in-plane spin component, which has helical spin texture for the states at 
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FIG. 5: (Color online) The schematic picture for the NaCo02 /.s-wave superconductor heterostruc- 
ture where Majorana bound state is associated with the vortex core. The chiral Majorana edge 
state is also expect at the edge of interface. 

Fermi level (see Fig. 4). The in-plane component is actually rather strong, and contribute to 
more than 10% of the total moment from our first-principles calculations. Considering the 
single dig state of NaCo0 2 surface and the 3-fold rotation symmetry, an effective surface 
Hamiltonian can be constructed as, 



(1) 



where £(k) = Ak + k_ + B(k+ + k Q _) + C(ArL&i) (with k± = k x ± iky) gives the non-spin- 
polarized band structure, V z a z is the exchange splitting, and the last term is the Rashba 



type SOC due to the surface. Evaluating the eigen values e(k) = £(k) ± JV* + a 2 k 2 with 
the parameters given in the caption of Fig. 4, the Fermi surface and its helical spin texture 
can be well reproduced (Fig.4(b)). Please note the hole type carrier and the negative sign 
of Rashba coupling a in our present case. 

When the NaCoC>2 surface comes into contact with an s-wave superconductor, a pairing 



term = ^k^c^c^ + h.c] will be generated by the proximity effect, and the 



ull Hamil- 



tonian reads H = H Q + which has been carefully studied previously |8Hl ll. Il9j . The 
dominant pairing channel should have the spin polarized p x +ip y symmetry, whose order pa- 
rameter is given as At(k) = -<*fc = k y +ik x ^ j n ^ limit of large V~, following Ref. [9fl, 

^ FK ' y/(V?+a 2 k 2 ) k to 25 to !_1J) 

the superconducting gap can be estimated as, 
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E„ 



(2) 



which is sizable. Here we use the parameters, effective mass m*=5mo, ol = — 0.066eVA, 
and V z =0.22eV, obtained from our first principles calculations. The size of A is determined 
by the superconducting gap of s-wave superconductor and its interface coupling (hopping) 
with NaCoC>2 side. 

Once the p x + ip y superconductivity is realized through proximity effect, the Majorana 
bound states associated with the vortex core will be expected [?], a, [l9|. The schematic 
experimental setting for the detection of Majorana fermions is shown in Fig. 5, where the 
absence of magnetic insulating layer and/or magnetic field is the essential difference com- 
pared to earlier proposals ffl. In additional, due to the broken time reversal symmetry 
(different from Fu and Kane's proposal |2j), a chiral Majorana edge state should exist at 
the interface (Fig.5). Since the single band half-metallic character of NaCo02 surface can 
be realized in a large energy window (of about 0.2eV), the fine tuning of chemical potential 
is not necessary. The large Fermi surface size and the high carrier density should be also 
helpful for sizable proximity effect, against the possible localization due to disorders. Among 
several possible choices of the superconductor, it is particularly interesting to consider the 
Na x Co02 superconductor with x around 0.35 [20|, [2jJ, since the lattice structures are well 
matched. The NaCoC^/s-wave superconductor heterostructure, therefore, provide us a new 
and simple platform for realizing topological superconductivity and Majorana fermion bound 
states. This work can be generalized in several directions. Similar predictions can be also 
made for AF Na^CoC^ (x >0.65) thin film with odd number of layers (which contribute to 
odd number of Fermi surfaces), or for LiCo02, which has the same crystal and electronic 
structure as NaCo02 [33]. 
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